C/(1)-Symmetry breaking and violation of axial symmetry in TlCuCla and other 
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We describe the Bose-Einstein condensate of magnetic bosonic quasiparticles in insulating spin 
systems using a phenomenological standard functional method for T = 0. We show that results 
that are already known from advanced computational techniques immediately follow. The inclusion 
of a perturbative anisotropy term that violates the axial symmetry allows us to remarkably well 
explain a number of experimental features of the dimerized spin-1/2 system TlCuCla. Based on 
an energetic argument we predict a general intrinsic instability of an axially symmetric magnetic 
condensate towards a violation of this symmetry, which leads to the spontaneous formation of an 
anisotropy gap in the energy spectrum above the critical field. We, therefore, expect that a true 
Goldstone mode in insulating spin systems, i.e., a strictly linear energy-dispersion relation down to 
arbitrarily small excitations energies, cannot be observed in any real material. 

PACS numbers: 75.10.Hk, 75.45.+j, 75.30.Gw, 73.43.Nq 



I. INTRODUCTION 

The concept of Bose-Einstein condensation (BEG), i.e., 
the occupation of a single quantum state by a macro- 
scopic number of bosons, has been extended from real 
bosonic particles to various types of quasiparticles with 
integer total spin [J, i, i, H i, i, H S i, 113, EH ■ 
Such discrete magnetic, electronic or lattice excitations 
are then treated as a set of bosons forming a Bosc gas 

. These quasiparticles usually possess a small effective 
mass, which permits to study BEG even at room temper- 
ature 11 1 . Both experiment [13, [l3| and theory [Tsl. [T3| 
suggest that the occurrence of a BEG in a 3D interacting 
Bose system has its origin in the spontaneous breaking 
of the U{1) symmetry, thereby leading to a gapless and 
linear excitation spectrum in the long- wavelength limit, 

1. e., to a Goldstone mode [l^. 

Experimental observations in a number of quantum spin 
systems can be explained within the theory of BEG, e.g., 
by a condensation of triplet states in dimerized spin- 
1/2 systems (hereafter called triplons) such as TlGuGla 
[1, 13, [l3l . The bosonic character of the these magnetic 
quasiparticles allows to describe this spin-dimer system 
as a weakly interacting Bose gas. Inelastic neutron- 
scattering measurements in the condensate phase of 
TlGuGla revealed, in accordance with theoretical inves- 
tigations [l3|, the presence of a seemingly gapless and 
linear excitation spectrum down to very low excitation 
energies of the order of 0.75 meV in ^o-f^ = 14 T [1]. 
This observation has been interpreted as a manifestation 
of the Goldstone mode. 

In the last years various scenarios for the consequences 
of anisotropy on the properties of a r nag netic BEG have 
been discussed in detail [H, [H, [13, [Sf [H . The pres- 
ence of any kind of anisotropy will, in principle, explicitly 
break the rotational (i.e, axial) symmetry of the bosonic 



system [2l|. The degree of spontaneous t/(l)-symmetry 
breaking then depends on the order of magnitude of 
the anisotropic terms compared to the energy scale as- 
sociated with the isotropic interactions of the Bose gas 
22 1 . The relatively large triplon bandwith in TlGuGla 
23| exceeds the spin-nonconserving terms, such as an in- 
tradimer exchange (IE) anisotropy and a Dzyaloshinsky- 
Moriya (DM) anisotropy, by more than two orders of 
magnitude [3, [H, [IH ■ Nonetheless, the question to what 
extent the existing anisotropics in TlGuGls do affect the 
magnetic phase diagram, the Goldstone mode and other 
measurable quantities of TlGuGla is still an issue under 
investigation [H, [H, [lU, [H . 

By taking a perturbative anisotropy term into account we 
will consider in the following the influence of an IE-like 
anisotropy that explicitly violates the axial symmetry, 
and we will study the consequences on the condensate 
phase of TlGuGk. The influence of a possible DM-like 
anisotropy [13, [IH is not considered here [1H|. Based 
on an energy consideration we will then argue that, as a 
consequence of an unavoidable magnetoelastic coupling, 
even an axially symmetric magnetic system is unstable 
towards a spontaneous violation of this symmetry as soon 
as the BEG-state is formed. 



II. FUNCTIONAL METHOD 

We describe the condensate at T = with a macro- 
scopic wave function, a complex scalar field ^{r, t). Stan- 
dard functional methods used to describe a dilute Bose 
gas in the classical limit at T = yield an extremal con- 
dition for the potential energy per dimcr (l^ , namely 



u{ip) — — H — - (^p^ip). = min. 

2 '11/1=1/11-1 
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where /i is the chemical potential, vq is a constant re- 
lated to a repulsive short-range interaction p^ . and tp^ 
is the complex conjugate of ip. The minimum value ipo 
then determines the condensate fraction ndO) = V'oV'Oj 
here defined as nc(0) = Nc/Nd with iVc the number of 
condensed triplons and Nd the number of dimers. 
In a dimerized antiferromagnet, the chemical potential 
is /i = gfiBfJ-aiH — He), where /is is the Bohr magne- 
ton, g the Lande g-factor, and He the critical magnetic 
field beyond which a triplet S ~ 1 state is energetically 
equally favorable as the singlet 5 = state. This can 
be expressed in terms of the energy gap A = g^g^oHc 
separating at zero field the S = 1 and the 5 = states, 
respectively. 

In the case of an explicitly violated axial symmetry we 
may include a perturbation term I7I (ipip + ip'^ip'^) jll. |2l| 
to the potential energy that can arise in a real magnetic 
system from various sources such as anisotropic intra- 
and interdimer interaction constants J and J, respec- 
tively For TlCuCla we have, for example, fioHc ~ 5.6 
T 1^, vo/kB = 315 K [13, and a I7I of the order of 
0.01 meV [H, [U, depending on the orientation of the 
magnetic field H with respect to the crystal lattice [l^ . 
Such an anisotropy term may arise from a pre-existing 
violated axial symmetry of the system, or from a spon- 
taneous distortion at the magnetic phase transition that 
we will discuss below. We, therefore, have to minimize 

u{ip) = -^i^H' + \i\ {i'yj + iJ^^^) + y {ijH'f , (2) 

where we assume that He itself is at first unchanged 
by the presence of the small perturbative anisotropy 
I7I < ^■ 

We first want to compare the results of this simple 
formalism with corresponding predictions from ad- 
vanced Hartree-Fock (HF) computations and with 
experimental data on TlCuCla. Despite the simple 
formalism used here, we can later make specific pre- 
dictions that would otherwise be more difficult to obtain. 



III. RESULTS 

A. Comparison with results from Hartree-Fock 
calculations 

Without any explicit anisotropy (i.e., \'y\ = 0) we ob- 
tain the well-known minimum value i/'oV'o = ndO) = /i/wg 
0. The phase of i/'o = |V'o|e"^ is not fixed in this case, 
leading to the u{ip) landscape sketched in Fig. [T] (left 
panel, "mexican-hat potential"). However, any nonzero 
value for I7I locks the phase of ipo to the imaginary axis 
(i.e., = ±7r/2 and therefore tpo = const.), and leads to 
an optimum value nc(0) = {fi + 2\j\)/vo in a "Napoleon's 
hat potential" (see Fig. [TJ right panel). The minimum 
potential energy becomes Umin = — (m + 2|7|)^/2uo which 
is smaller than in the axially symmteric case. The 




FIG. 1: Potential-energy it as a function of ip for axial sym- 
metry (left panel, "mexican-hat potential") and violated axial 
symmetry (right panel, "Napoleon's hat potential"), respec- 
tively. In the symmetric case the minimum value Umin ~ 
—fi^ /2vo is realized along a circle (dashed line), while in the 
anisotropic case isolated minima Umi»i = — (p-|-2|7|)^/2wo are 
on the imaginary axis (filled circles). 



saddle-point value for u on the real axis for fi > 2|7| 
is ~{fi — 2|7|)^/2wo- Note that these energy densities are 
all expressed per dimer. The corresponding energies per 
triplon are u/nc(0). 

A vanishing ndO) is realized when {fi + 2|7|) = 0. As 
a consequence, the gap field He that would be observed 
in an ideal system with axial symmetry is renormalized 
to a value H^^'p = H^- AH^ (with AH^ = 2\j\/ gfiofJ-s) 
above which condensation occurs. Takin g a reasonable 
value for I7I w 0.01 meV for TlCuClg and H || b 

with g =2.06 [l^, we obtain a renormalization of the 
critical field due to I7I alone by ^qAHc « 0.2 T for this 
particlular magnetic-field direction. 

The resulting condensate fraction, ndO) = fi'^^^/vo where 
^eip _ gfiBfJ,o{H — H^^P), is in full agreement with 
Hartree-Fock calculations for spin dimer systems In 
Fig. [2^ we show the triplon condensate fraction ndO) 
as deduced from our magnetization M(T, _ff)-data of 
TlCuCls [11], see Fig. that we have already cor- 
rected for a small fraction h of noncondensed triplons 
[^ . These data have been obtained from the simple rela- 
tion M(0) = g//Bn(0)A'(i (with Ad the number of dimers 
and n(0) the total triplon fraction at T = 0) with- 
out assuming any specific value for vq . A linear fit to the 
data in the dilute hmit {^qHc < f^oH < 7.5 T [U) yields 
Aioi^r^ = 5.501± 0.003 T and vo/kB = 311.4 ± 0.5Km3 
which is in very good agreement with available litera- 
ture values [111, [13, [1^. The deviation from the linear 
behavior at larger magnetic fields can be attributed to 
the contribution of higher triplet states [^. Our simple 
formalism does not include the infiuence of such triplet 
states, nor does it allow for a determination of h itself, 
but this latter correction is of the order of a few percents 
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FIG. 3: a) Potential energy per dimer, and b) energy-level 
scheme of the lowest triplon branch of an axially symmet- 
ric system (I7I — 0, thick solid lines), and of a system 
that shows a spontaneous axial distortion (I7I > 0) above 
H* (thin lines and arrows, axes are not to scale), respec- 
tively. The grey shaded area indicates the maximum hystere- 
sis width of a possible weakly first-order transition at H* with 
jjexp ^ j^* ^ jj^ (filled arrows: H increasing; open arrows: 
H decreasing) . The dashed line in b) reproduces the trend in 
the ESR data from Refs. [3 and [Til for TlCuCla and H || 6 
(see text). 



FIG. 2: a) Condensate fraction nc(0), non-condensed triplon 
fraction n(0), and the percentage of non-condensed triplons 
(inset) [2^ . b) Triplon-fraction n(0) data obtained from mag- 
netization M(T)data for H || fe (see arrow for the fioH = 9 T 
data and Ref. [2g]). c) Schematic representations of the ef- 
fects of violated axial symmetry on the triplon fraction at 
T = and at T = T*^^, d) on the triplon fraction at fixed 
magnetic field H > He [1, [20l . [2ll | and e) on the excitation 
spectrum E{k) (axes are not to scale). Solid and dashed lines 
represent an axially symmetric system and a system with vi- 
olated axial symmetry, respectively. 



at most in our data [3|, |28[ , as it is typical for a weakly 
interacting Bose gas, see inset of Fig. [5^. 
Focusing further on the effects of an explicit violation 
of axial symmetry, the condensate fraction at T = 
is changed by 2171/110 in our calculation, see Figs. 2c 
and d. If we use the fact that the total triplon fraction 
n(0) « nc(0) and take the result n{T^^) = n(0)/2 from 
Ref. 3 assuming a quadratic triplon dispersion relation, 
we obtain a shift |7|/t'o in ^{T^^) that is again in full 
agreement with the corresponding HF calculations [2]| . 



see again Figs. 2c and d. 

We may relate the minimum value Umin at T = 
to the transition temperature T^^ if we assume that 
\umin\/nc{Q), the energy gain per triplon upon conden- 
sation, is proportional to ksT^^ with a field-dependent 
proportionality factor of the order of unity. Any nonzero 
I7I leads to an increase in T^^ as compared to the axially 
symmetric case, see Fig. This trend can be cle arly 
seen in the calculated Af (T, iJ )-curvcs from Ref. [2l| , 
where the minimum in AI that is usually taken as a crite- 
rion to define T^^ is shifted towards higher temperatures 
as soon as I7I > 0. If we again take I7I = 0.01 meV we ob- 
tain a shift in T^^ of the order of AT^^^ « \l\/kB « + 
O.IK, which has to be compared to the result of the more 
precise HF calculations with AT^^^ « + 0.5 K [H. 
As a consequence of the anisotropy term I7I the orig- 
inal invariance of u(V') with respect to a transforma- 
tion ■!/; — > ijje^'^ is lifted. The maximum variation in 
the potential-energy per dimer u along the ellipsoid con- 
tour with local minima in the radial |'!/'|-direction, see 
Fig. [U is 4/x|7|/wo, or A w 4|7| per ground-state triplon. 
This determines the order-of-magnitude of an anisotropy 
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gap that can be calculated to ^ = ^8|7|g^_B/io(i? — He) 

[Hi, [131. 

Excitations above the ground state tpQ with energies be- 
low A will clearly not show the typical gapless Goldstone- 
like behavior as expected for an axially symmetric sys- 
tem. For TlCuCla we calculate with I7I = 0.01 mcV 
[li,[2l| a gap z4 w 0.3 meV for ^i^H = 14 T along 6, which 
is somewhat below what has been resolved in inelastic- 
neutron scattering measurements [3|. For excitation en- 
ergies larger than A the presence of an anisotropy gap 
may remain unnoticed, see Fig. [5^. 

In TlCuCls such a gap may arise from a pre-existing 
anisotropy that is already present in = [l^. In 
the following we argue, however, that even a perfectly 
axially symmetric magnetic system is unstable towards 
a spontaneous violation of this symmetry at H^, which 
inevitably leads to the formation of a small anisotropy 
gap A above He of real materials. 



B. Instability of the condensate towards violation 
of axial symmetry 

A striking fact in our analysis is that the minimum 
potential energy per dimer is smaller with a nonzero I7I 
than in an analogous axially symmetric system with I7I 
= 0. This means that a distortion of the original crystal 
symmetry may spontaneously occur at He together with 
an increase of I7I, provided that the total energy, includ- 
ing both magnetic and crystal-lattice contributions, is 
lowered along with this distortion. This argument is so 
general that it should be applicable to all insulating spin 
systems that are supposed to show a Bose-Einstein con- 
densation of magnetic bosonic quasiparticles. As we do 
not not make any specific assumptions on the microscopic 
arrangement of the spin-carrying atoms, one cannot make 
any more precise, universal statement about the details 
of the resulting lattice distortion. 

The gain in potential energy per dimer upon condensa- 
tion in combination with this simultaneous distortion is 
2|7p/wo, see Fig. [3^. If the critical field is approached 
from below with H increasing, the parameter I7I may 
therefore jump discontinuously from zero to its opti- 
mum value either at He or at a transition field iJ* with 
Ijexp £ ^* £ jj^^ while with H decreasing the transition 
can take place at a different field in the same magnetic- 
field interval, see grey shaded area in Figs. [3^ and b. 
In an ideal situation with a perfect axial symmetry in 
H = 0, the critical field He corresponds to a "normal- 
state" value with I7I = that is determined only by 
the gap energy, while in the condensate phase the effec- 
tive critical field is H^^p = He — AHe < He, lowered 
by AHe = 2|7|/g//sMo with respect to He due to the 
increase in |7|. If the transition occurs at a transition 
field H* that is strictly larger than H^^^ , one will ob- 
serve at T = corresponding small discontinuities in u 
(Au < 2|7|2/wo, see Fig. W), (An^ < 2|7|/wo, see 
Fig. [!:) and M (AM = AuegfisNd), which would qual- 



ify the transition as of weakly first order with a maxi- 
mum observable hysteresis width AHe- The occurrence 
of hysteretic effects in a real material may depend, how- 
ever, on further conditions that are not considered here, 
such as material-quality issues or the relevance of possi- 
ble quantum fluctuations at T = to the order of the 
phase transition. 

Even if a non-zero I7I evolves continuously above the crit- 
ical field as a function of H along with a continuous 
structural distortion with no detectable hysteresis, one 
should still be able to distinguish between the normal- 
state He extracted from the experimental data taken be- 
low He^P, and a He^^ < He from corresponding data 
taken well above He, respectively, to obtain an estimate 
for I7I from the resulting difference AHe- In any case, 
if He^P < H* < He, one expects a finite energy dif- 
ference < I.6I7I at the transition field H* between the 
lowest-triplon state and the singlet states, respectively, 
see Fig. [SJd. 

A possible pre-existing anisotropy that may already be 
present at if = (which is likely to be the case in 
TlCuCls) can be easily included in this formalism by 
identifying He with a renormalized " normal-state" value 
that already contains this pre-existing anisotropy. Any 
additional I7I that may develop together with the lattice 
distortion around this critical field will somewhat change 
the value of the anisotropy gap A. The difference AHe, 
however, and possible discontinuities in u, Ue and M are 
determined by the additional I7I alone. 
It is remarkable that electron-spin resonance (ESR) data 
taken on TlCuCla do indeed show a clearly gapped be- 
havior at and above H = He^P \\ b [l8l.[l9j. as we sketched 
in Fig. [3)3. Moreover, the ESR frequencies due to the low- 
est triplon gap in the "normal" phase for H || 6 extrapo- 
late to zero at a somewhat larger He (by fioHe ~ 0.2 T) 
than the square-root like gap that we attribute to A{H) 
in the condensate phase and that can be fitted nicely with 
a I7I « 0.016 meV, see energy-level scheme in Fig. [3)d. 
Our scenario may also be an explanation for the observed 
abrupt changes in the "^^Cl quadrupole shift that has 
been interpreted as an indication of a weakly-first order 
lattice deformation, as well as for the pronounced hys- 
teretic behavior (with a fi^AHe ~ 0.2 — 0.3 T) of the 
observed peaks in the sound-attenuation data of TlCuCla 
at Heir) [3l|. These observations may indicate that I7I 
of TlCuCla is indeed larger for H > He than well below 
this value. 

The present picture may also account for the first-order 
like features that have been seen in the X-ray data of the 
spin-ladder compound Cu2(C5Hi2N2)2Cl4 [S^] and in the 
magnetocaloric effect of the axially symmetric S = 1 sys- 
tem NiCl2 — 4SC(NH2)2 [13 at the respective magnetic 
phase transitions. It is also not unreasonable to assume 
that the observed gap feature in the ESR-data of this 
latter compound at fioH = 8 T [s^l is also related to a 
possible lattice distortion. 

Such a spontaneous distortion arising from the interplay 
between the emerging magnetic Bose-Einstein conden- 
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sate and its host crystal that lowers the total energy, with 
a tendency to increase or even create an anisotropy per- 
pendicular to the external magnetic field above even 
in a perfectly axially symmetric system, is rather unique 
and is reminiscent of the spin-Peierls instability in one di- 
mensional magnetic chains. A spin-Peierls-type scenario 
has indeed been suggested to explain the NMR data of 
the spin-ladder compound Cu2(C2Hi2N2)2Cl4 around its 
critical magnetic field [H, [s^. This instability should 
be a universal feature of magnetic BEC systems at their 
magnetic phase transition, and it is not expected to occur 
in axially symmetric Bose gases composed of real parti- 
cles such as superfluid "'He or atomic condensates, where 
the condensate cannot create an axial-symmetry break- 
ing term by itself. 

In a microscopic picture, the tendency of a magnetic con- 
densate to spontaneously violate the axial symmetry can 
be interpreted as a natural consequence of the trans- 
verse magnetic ordering that develops in the condensate 
phase and that locks to the crystal lattice due to un- 
avoidable magnetoelastic coupling [3, [13, IH]. As soon 
as the transverse magnetic moments point to a specific, 
energetically preferred crystal direction, the phase that 
is associated with the angle between these moments and 
the crystal axes [iBl is indeed fixed (in TlCuCla with an 
angle a « 39° to the a-axis [131), ^'^'^ ^^e magnetic ana- 
logue to a supercurrent velocity = h/m*V(j> (where 
m* is the effective mass of a triplon) is naturally zero for 
excitation energies below the anisotropic gap A. This 
gap covers an excitation-frequency range that is crucial 
for experiments that rely on the existence of a long-lived 
phase coherent condensate [^, such as the detection of 
a long-lived spin supercurrent as in '^He — B [1^ of 
macroscopic second-sound like oscillations as observed in 
superfiuid ^He [4l| or of stable vortex-like structures as 
they have been observed both in superfluid ^He and in 



atomic condensates [42|. To achieve a lifetime of the 
order of seconds for a phase coherent condensate, a cor- 
responding anisotropy gap may not exceed a few fcmto 
electron volts. 

In the absence of a Goldstone mode (for which all values 
of the phase have to be energeticaly equivalent), the 
quantity "00^0 in the zero-frequency limit docs no longer 
represent a condensate fraction, but it is rather related to 
the order parameter characterizing the antiferromagnetic 
state [l|. 



IV. CONCLUSIONS 

We have analyzed the spontaneous symmetry break- 
ing in a Bose gas of magnetic bosonic quasiparticles 
in insulating spin systems based on simple functional 
methods in the classical limit. Our results reproduce 
several results from earlier HF approximated computa- 
tions [3, [20l . [2ll |. and various experimental findings in 
TlCuCla such as the occurrence and the magnitude of 
an anisotropy gap [H, [3 a.nd a weakly first-order like 
behavior at the magnetic phase transition [U, [3l| can 
be explained. On the basis of an energetic argument we 
expect that all magnetic BEC systems in insulating spin 
systems are intrinsically unstable towards a spontaneous 
anisotropic distortion perpendicular to the external mag- 
netic field, which leads to the formation of an anisotropy 
gap that is seriously limiting the lifetime of a phase co- 
herent condensate. 

We thank to O. Sushkov, A. K. Kolezhuk and E. Ya. 
Sherman for stimulating discussion. This work was 
partly supported by the Schweizerische Nationalfonds zur 
Forderung der wisscnschaftlichcn Forschung, Grant. No. 
20-111653. 
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